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The lunar Askaryan technique, which involves searching for Askaryan radio pulses from particle 
cascades in the outer layers of the Moon, is a method for using the lunar surface as an extremely 
large detector of ultra-high-energy particles. The high time resolution required to detect these 
pulses, which have a duration of around a nanosecond, puts this technique in a regime quite 
different from other forms of radio astronomy, with a unique set of associated technical chal¬ 
lenges which have been addressed in a series of experiments by various groups. Implementing 
the methods and techniques developed by these groups for detecting lunar Askaryan pulses will 
be important for a future experiment with the Square Kilometre Array (SKA), which is expected 
to have sufficient sensitivity to allow the first positive detection using this technique. 

Key issues include correction for ionospheric dispersion, beamforming, efficient triggering, and 
the exclusion of spurious events from radio-frequency interference. We review the progress in 
each of these areas, and consider the further progress expected for future application with the 
SKA. 
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1. Introduction 

The lunar Askaryan technique, first proposed by Dagkesamanskii & Zheleznykh [2], employs 
the Moon as a detector of ultra-high-energy (UHE) particles by searching for the Askaryan radio 
pulse [1] produced when such a particle interacts on the lunar surface. The large size of the Moon 
gives such a detector a very large aperture, but the distance to the Moon makes the expected pulse 
extremely weak, requiring a highly sensitive radio telescope to detect it. The Square Kilometre 
Array (SKA) [3] will be, when complete, more sensitive in the appropriate frequency range than 
any current radio telescope by around an order of magnitude [4], making it highly attractive for this 
application. This — and the precision measurement of air showers with the same instrument [5] — 
is the focus of the SKA High Energy Cosmic Particles (HECP) Eocus Group'. 

Past experiments [6, 7, 8, 9, 10, 11, 12, 13, 14] have been sensitive only to the UHE neutrino 
fluxes predicted by certain exotic ‘top-down’ UHE particle production models (e.g. [15, 16]). While 
the SKA will not be sufficiently sensitive to detect the lower-energy cosmogenic neutrino flux, it 
will be able to further constrain top-down models — and, critically, it will be able to detect the 
known UHE cosmic ray (CR) flux, providing the first positive detection of a UHE particle with 
this technique, and allowing it to test the isotropy of the UHECR sky. Eor further details of the 
scientific potential, see the co-contribution by James et al. [17]. 

In this contribution, we discuss some of the issues faced by this experiment, the experience 
gained from past efforts, and the further work to be done before they are overcome. Section 2, 
contains a general description of the relevant chai'acteristics of the SKA, while section 3 contains 
more details of its beamforming capabilities. In section 4, we describe the capabilities required to 
correct for dispersion of an Askaryan pulse as it passes through the ionosphere, and in section 5 
we address the issue of spurious signals from radio-frequency interference. Einally, section 6 sum¬ 
marises the outlook for overcoming these issues for successful application of the lunar Askaryan 
technique with the SKA. 

2. The Square Kilometre Array 

Phase 1 of the SKA, to be constructed in 2018-2023, will consist of two components. SKAl- 
EOW will be an electronically-steered aperture array, similar to the pathfinder instruments EO- 
EAR [19] and MWA [20], with a frequency range of 50-350 MHz and situated in Western Aus¬ 
tralia. SKAl-MID will be an array of steerable parabolic dishes, with a selection of receivers to 
cover the frequency range 350 MHz to 13.8 GHz, located in South Africa. Phase 2, to be con¬ 
structed in 2023-2030, is not yet so well defined, but will most likely consist of a large expansion 
of both of these instruments. 

Past lunar Askaryan experiments have used a range of frequencies, from 113 MHz [12] to 
4.8 GHz [8], trading off between the greater Askaryan pulse strength (and hence lower minimum 
detectable particle energy) at high frequencies and the increased opening angle of the Askaryan 
emission (and hence greater geometric particle aperture) at low frequencies [21, 22]. The HECP 
Eocus Group is focusing exclusively on SKAl-EOW, taking advantage of the latter effect. 

'http://astronomers.skatelescope.org/home/focus-groups/high-energy-cosmic-particles/ 
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Figure 1: Beamforming strategy for lunar Askaryan observations with the SKA, with the beam from each 
station centred on the Moon, and all core stations synthesised into multiple array beams directed at different 
points on the Moon. These points are selected around the limb of the Moon, from which a signal is most 
likely to be detected [24]. The size of the station beam shown is for the FWHM at the top of the 100- 
350 MHz band, and the array beams are for the geometric centre of this range. 

SKAl-LOW will consist of 131,072 log-periodic dipole antennas [23], grouped into ~ 450 
stations of 35 m diameter, each with ~ 300 densely-packed antennas. The majority of the stations 
will be grouped within 1 km, forming a dense core to the array, with the remainder spread out over 
baselines of up to ~ 80 km. The total collecting area will be less than 1 km^; this milestone will be 
reached with phase 2. 

3. Beamforming 

Because the antennas of SKAl-LOW are immobile, pointing of the telescope will be done 
electronically, with signals from all antennas in a station being summed with appropriate delay 
offsets to form a station beam pointed at an arbitrary point on the sky. The full width half maximum 
(FWHM) size of a station beam has a minimum value of 1.4° at the top of the band, which is 
sufficient to view the entire Moon, as shown in figure 1 . 

For high-time-resolution observations such as these, SKAl-LOW will have an array-level 
beamformer, combining multiple station beams in a similar fashion to form an array beam. By 
summing station beams with different delay offsets, this beamformer will be able to form multiple 
array beams simultaneously, with a planned capacity of 16 dual-polarisation beams with the full 
bandwidth. The size and sensitivity of these array beams depends on the set of stations used to form 
them. If only stations in the dense core of the array are used, the beams are large enough for them 
to be tiled completely around the limb of the Moon (see figure 1), achieving an efficiency facfor of 
~ 50% of fhe particle aperture for full lunar coverage, which is incorporafed info fhe simulations 
of [17]. The sensitivity is reduced compared to the entire array, but if the array beams are used to 
trigger storage of buffered data from all stations (per [18]; see figure 2), fhe full sensifivify of fhe 
array can be achieved refroacfively. 
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Figure 2: Preliminary design of the signal path for a lunar Askaryan experiment with the SKA. The entire 
signal path (upper left) includes transient buffers, to be stored when triggered by the detection of a candidate 
pulse. The triggering unit (lower left) performs the real-time search for candidates, with each beam being 
processed by a search module (right). 


The application of a beamformer to the lunar Askaryan technique has already been demon¬ 
strated in the NuMoon WSRT experiment [12]. However, the case of SKAl-LOW differs in that 
the inputs and outputs of this instrument’s beamformer are in the form of channelised data produced 
by a polyphase filter (PPF), so the PPF channelisation must be inverted to restore the original high- 
time-resolution signal. This inversion has been shown to be possible and adequately efficient as 
part of preparatory work for a proposed future NuMoon experiment with LOFAR [25]. 

Retroactively achieving the full array sensitivity requires beamforming in software with stored 
data. This capability has been demonstrated, for the detection of slower astronomical transients, 
by the FRATs project [26]. Of lunar Askaryan experiments, the only one to attempt this thus far 
is the LUNASKA Parkes-ATCA experiment, which completed some of the required calibration 
steps [14], but has not demonstrated the capability to form retrospective beams and search them for 
Askaryan pulses. 
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Figure 3: Left: a radio pulse (solid) with a flat spectrum across the frequency range 100-350 MHz, and the 
same pulse after dispersion equivalent to an electron column density of 0.1 TECU (dashed) and 1.0 TECU 
(dotted). Right: the fraction of the original pulse amplitude recovered after this dispersion (or dedispersion 
error), for several frequency ranges. A spacing of ATEC = 0.1 TECU between dedispersion trials ensures 
that the true dispersion will lie no further than 0.05 TECU from the closest trial, giving a typical recovered 
pulse amplitude of > 95% for any of the frequency ranges shown here. 


4. Ionospheric dispersion 

A lunar-origin radio pulse propagating through the ionosphere experiences a frequency-dependent 
delay. This ionospheric dispersion extends the pulse in time and decreases the peak pulse height, 
as shown in figure 3 (left). The degree of dispersion is determined by the column density of free 
electrons, or total electron content (TEC), measured in TEC units (TECU; 1 TECU = 10^® elec¬ 
trons m^^). It also depends on the frequency range of the pulse: the low frequency and broad band 
of SKAl-EOW make it particularly severe. 

Previous experiments have either neglected dispersion, or corrected for the dispersion associ¬ 
ated with the measured ionospheric TEC. Eor SKAl-EOW, this is insufficient: since the recovered 
amplitude of the pulse drops so quickly with any error in the TEC (see figure 3 (righf)), even 
affer excluding fhe 50-100 MHz band which is most strongly affected by dispersion, a typical un¬ 
certainty in the measured TEC of 1-2 TECU [13] can result in a loss of over 50% of the signal 
amplitude. 

A direct solution to this problem is to search in parallel at multiple TECs, as shown in figure 4 
(left), at the expense of increased processing requirements. To span an uncertainty range of ±1- 
2 TECU with a spacing between filters of ATEC ~ 0.1 TECU would require ~ 30 filters running in 
parallel. If implemented as finite impulse response (EIR) filters, each of these would require ~ 100 
taps (operations per sample), per figure 4 (righf). This mofivafes fhe 32 parallel EIR fillers wifh 
128 taps each in the preliminary search module design shown in figure 2 (right), which dominate 
the total computing requirements. 

The computational requirements can be considerably reduced if the TEC can be measured 
more precisely: the number of parallel filters, and the number of operations per filter, are both 
decreased. The quoted uncertainty of 1-2 TECU is for two-dimensional TEC maps based on data 
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Figure 4: Left; planned dedispersion scheme, with bulk dedispersion compensating for the mean expected 
TEC, and fine dedispersion to implement a parallel search with spacing ATEC across the uncertainty range 
defined by ciTEC. Right; required length of an EIR dedispersion filter to retain > 98% of the power of a 
dispersed template pulse. To accommodate a dispersion equal to the uncertainty of crTEC ~ 1 TECU, the 
required filter length is ~ 100 taps, depending on the frequency range. 


from global positioning system satellites; direct measurements of the TEC along the line of sight 
to these satellites have an inherent uncertainty of <0.1 TECU [27], and direct interpolation of 
these measurements should allow something closer to this to be achieved. Another possibility is 
to use Earaday rotation of polarised lunar thermal radio emission, which is also determined by the 
ionospheric TEC, to directly measure the TEC along the line of sight from the telescope to the 
Moon [28]. 


5. Radio-frequency interference 

A critical problem for a lunar Askaryan experiment is to apply cuts to exclude spurious an¬ 
thropogenic radio pulses, referred to as radio-frequency interference (REI), which can appear as a 
false signal. In application to the SKA, this problem comes in two parts. The first is to exclude 
in real time the majority of REI, which would otherwise cause an excessively high trigger rate 
and an unmanageable quantity of stored data. The second is to completely exclude any remaining 
REI in retrospective analysis, so a detected pulse can be identified wifh high confidence as being a 
lunar-origin Askaryan pulse produced by a UHE particle. 

Pasf experimenfs have mosf commonly used a coincidence requiremenf falling info one of fwo 
cafegories; requiring a coincidenf defecfion by mulfiple defecfors poinfed af fhe same parf of fhe 
Moon, or requiring an anficoincidenf defecfion by only one of multiple defecfors poinfed af differenf 
parfs of fhe Moon. The former approach is inapplicable here, so we are forced fo rely on fhe second, 
wifh fhe 16 array beams consfifufing our multiple defecfors. The pasf experimenfs fhaf have relied 
mosf heavily on fhis approach (alongside ofher cufs) are fhe NuMoon WSRT experimenf [12] and 
fhe EUNASKA Parkes experimenf [13]. 

In fhe NuMoon WSRT experimenf, a subsfanfial number of REI evenfs passed fhe cufs, domi¬ 
nating sfrongly over fhe fhermal noise disfribufion. This may be due fo fhe REI-loud environment 
of the telescope, the sensitivity of the array beams to the REI-loud horizon, or the low observing 
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frequency, which typically results in more RFI. The LUNASKA Parkes experiment had a higher 
observing frequency and a more RFI-quiet environment, and the cuts applied were sufficient to 
remove effectively all RFI, but the statistical rigour of this was not clearly established: if a single 
event had remained after the cuts were applied, it would not be clear if it was a remnant RFI event 
or a true lunar-origin pulse. 

Despite its low frequency, SKAl-LOW has two major advantages over these two experiments. 
The first is the telescope site, which is extremely isolated and therefore RFI-quiet. The second is 
the availability of station-level buffers which permit retrospective beamforming, through which it 
is possible to precisely localise a source of RFI in the near field of fhe telescope. Togefher, fhese ad- 
vanfages make if highly likely fhaf if will be possible fo meef fhe sfandard of proof required fo make 
a definitive detection of a lunar-origin pulse, buf esfablishing Ibis wifh cerfainfy will require analy¬ 
sis of stored dafa from SKAl-LOW in operafion. Until Ibis is possible, a reasonable intermediate 
sfep would be fhe measuremenf and analysis of fhe nanosecond-scaled pulsed RFI environmenf af 
fhe telescope sife. 

6. Conclusion 

The Square Kilomefre Array provides a pofenf option for applying fhe lunar Askaryan tech¬ 
nique, wifh a sfrong prospecf of ifs firsl posifive detection of an ulfra-high-energy cosmic parficle, 
buf furlher fechnical developmenf is necessary fo demonsfrafe fhe capabilify fo fully exploif Ibis in- 
sfrumenf. Experimenfs of Ibis fype have nol previously searched for Askaryan pulses in refrospec- 
fively beamformed dafa; fhis capabilify could be demonsfrafed fhrough a complefe analysis of fhe 
LUNASKA Parkes-ATCA experimenf, or fhrough fhe realisation of fhe planned NuMoon LOFAR 
experimenf. Currenf measuremenfs of fhe ionospheric elecfron confenf are sufficienl for dedisper¬ 
sion wifh fhe SKA, buf fhe compufing requiremenfs could be subsfanlially decreased fhrough eifher 
direcf inferpolafion of safellife measuremenfs, or measuremenf of lunar Faraday rofafion. The ef- 
fecliveness wifh which false signals from RFI can be excluded cannof be fully esfablished until 
an experimenf is carried ouf, buf exfrapolafion from previous experimenfs suggesfs fhaf if will be 
pracfical, and an on-sife survey of fasl fransienf RFI would help fo demonsfrafe fhis. 

None of fhese fechnical issues appears likely fo presenf an insurmounfable obsfacle fo lunar 
Askaryan observafions wifh fhe SKA. The fask now af hand is fo demonsfrafe solufions fo fhem 
before fhe scheduled complefion of fhe SKAl-LOW insfrumenf in 2023. 
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